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PREFACE 

In  this  RAND  Memorandum,  the  author  presents 
mathematical  results  on  radiative  transfer  in  slab  geometry. 
This  subject  has  important  implications  for  meteorology, 
astrophysics,  and  the  detection  of  nuclear  blasts. 

The  research  presented  here  was  sponsored  by  the 
Advanced  Research  Projects  Agency. 


-V- 


SUMMARY 

In  this  Memorandum  a  new  probabilistic  approach  to 
radiative  transfer  problems  is  presented  in  such  a  way 
that  the  integral  equations  for  the  stationary  and  non¬ 
stationary  scattering  functions  in  a  finite  inhomogeneous 
flat  layer  are  derived  directly  from  a  somewhat  modified 
form  of  the  Chapman-Kolmogorof f  equation. 
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A  PROBABILISTIC  APPROACH  FOR  SCATTERING  OF  LIGHT 
IN  SLAB  GEOMETRY. I. 


1 .  INTRODUCTION 

While  the  exact  theory  of  radiative  transfer  in  the 
initial  stage  has  been  mainly  developed  by  several  authors 
(Hopf  [8],  Ambarzumian  [l],  and  Chandrasekhar  [7]),  the 
following  powerful  methods  stemmed  recently  from  the 
principle-of-invariance  method:  the  invariant-imbedding 
technique  due  to  Bellman  and  Kalaba  [2]  (see  also  Wing  [ ll] ) , 
the  combined-operations  method  of  Busbridge  [5],  and  the 
probabilistic  method  of  Sobolev  [9]  and  Ueno  [ 10] . 

Of  the  above,  the  combined-operations  method  and  the 
probabilistic  method  do  not  permit  us  to  derive  directly 
an  integral  equation  for  the  scattering  function  in  slab 
geometry,  whereas  the  invariant-imbedding  technique  does 
so.  In  the  present  paper,  a  new  probabilistic  approach 
for  the  derivation  of  the  integral  equation  governing  the 
scattering  function  in  flat  layer  is  elucidated.  Even 
when  there  is  some  distribution  of  emitting  source  in  the 
medium,  we  can  compute  the  angular  distribution  of  emergent 
radiation  from  the  medium  by  means  of  the  scattering  function. 

The  characteristic  feature  of  the  present  approach  is 
that  an  integral  equation  for  the  scattering  function 
defined  as  a  truncated  generating  function  of  the  emission 
probability  is  derived  directly  from  the  Chapman-Kolmogorof f 
equation  in  somewhat  modified  form.  In  Sec.  2,  we  shall 
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treat  the  stationary  intensity  diffusely  reflected  by  a 
finite  inhomogeneous  flat  layer,  and  in  Sec.  3  the  time- 
dependent  diffusely  reflected  intensity  by  a  finite 
inhomogeneous  flat  layer  will  be  considered. 

2.  THE  STATIONARY  CASE 

Consider  a  plane-parallel,  inhomogeneous,  and 
isotropically  scattering  atmosphere  of  finite  optical 
thickness  7^  >  0,  when  the  lower  surface  r  ■  is 

illuminated  by  axially  symmetric  radiation  of  intrinsic 
flux  F  at  an  angle  cos-1fio(0  <  <  D  with  the  inward 

normal  to  the  surface  T  »«  . 

Let  I(T,  +  M)(0  <  4  <  1)  and  I  (T,  -  M)  (0  <  M  <  D 
denote,  respectively,  the  intensities  of  radiation  directed 
toward  the  surfaces  T  *  and  T  =  0. 

The  equation  of  transfer.  The  equation  of  transfer 
appropriate  to  the  present  case  takes  the  form 

(2.1)  =  Hr,  H)  -  /|l(r, 

where  o>  is  the  albedo  for  single  scattering,  i.e.,  the 
probability  of  photon  survival,  together  with  the 
boundary  conditions 


(2.2)  I  (0 ,  +  M)  =  0  and  I  (r^ ,  -  M)  =  tj  8  (M  -  IJ-Q)  (0  <  fi  <  1)  . 


NvmMMi  mkssmhmi*** 


Then  the  equation  of  transfer  in  the  diffuse  radiation 


field  is  expressed  in  the  form 


(2.3)  M 


dl(r, 


dr 


=  I  (r, 


M) 


£D(T) 

"5 


r+l 

J  Ur, 

-1 


jLt' ) d/x'  -  |a)(r)e_(Tl  /M 


Equation  (2.3)  should  be  solved  subject  to  the  boundary 
conditions 


(2.4)  1(0,  +  M)  =0  and  1(7^,  -  11)  *  0  (0  <  M  <  D  • 


The  emission  probability  distribution.  Let 
p*(M;  t,  TjMO  <  M  <  1»  0  <  T  <  T!>  denote  the  probability 
that  a  photon  absorbed  at  the  level  T  will  reappear  in 
the  direction  M  in  the  radiation  emerging  from  the 
surface  t  »  r^. 

The  integral  equation  for  p*(M;  T ,  r^)  is  written 
[10]  in  the  form 

(2.5)  p*(M;  T,  7^)  =  o)(r)[e_^Tl  ^  +  7VT{p*(Ml  t,  7^)}], 

where  the  truncated  Hopf  operator  h  is 

(2.6)  \{f(t)}  =  \  ^1f(t)E1(lt  -  Tl)dt. 

In  (2.6),  Ex(t)  is  the  first  exponential  integral 


defined  by 


JKimum  VMttW  *«T**«*W 
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,  .  r1  -t/4  dM 

(2.7)  E, (t)  =  J  e  u~  • 
x  0 

In  this  context  it  is  mentioned  that  a  photon 
diffusion  equation  (2.5)  is  derived  [10]  from  the  transfer 

equation  (2.3)  . 

The  scattering  function.  Let  the  truncated  generating 
function  of  the  emission  probability  distribution 
p *(M;  T,  T±)  be  denoted  by 

(2.8)  n (M/M0,  T!>  “  ^  p*(4:  T*  Tl)e  (Tl  T)/M°  dT‘ 

The  above  function  may  be  considered  also  a3  the  truncated 
characteristic  function  of  the  emission  probability  defined 
in  terms  of  the  bilateral  Laplace  transform. 

On  the  assumption  that  the  truncated  generating 
function  IXM  1  Tj.)  represents  also  the  Markovian 

conditional  probability  of  diffuse  reflection  at  the  lower 
surface  T  «  T1#  an  integral  equation  for  UQJ.  1  V 

follows  directly  from  the  Chapman-Kolmogorof f  equation  in 
somewhat  modified  form  given  by 

(2.9)  Ii(M  l  Ti  +  ^  ~  ^  *  ^o’  Tl^p^  *  :  Tl»  > 

for  0  <  A  <  Tx ,  where 

(2.10)  !T(M  l  md,  tx)  *  n(M  I  Tx)(l  -  jj~  )  +  A-  p*(M;  T,  tx), 
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and 


(2.11)  p(m  I  M’ ;  V  a)  =  2£,p*(m;  V  V  +  5(4  4’)d  -  4.) 


In  (2.10),  the  first  term  on  the  right-hand  side  is 
due  to  the  probability  of  diffuse  reflection  of  photon  by 
the  additional  layer  A  without  scattering,  and  the  second 
term  represents  the  probability  of  photon  emergence  from 
the  medium  with  scattering  in  the  additional  layer  A. 
Equation  (2.11)  is  equivalent  [ 10]  to  the  conditional 
transition  probability  used  in  the  derivation  of  an  integral 

equation  for  p*(M;  T>  T^) • 

Letting  A  -►  0  in  (2.9),  we  have 


(2.12) 


an 


+  ( 


4 


r)n 

o 


p*(4;  t 


ri)U 


5  /» 


(4’ 


T  )*£l' 

1  V 


From  (2.5)  we  get 


(2.13)  p *(4 


;  r1,  rx)  =  o)(r1)[l  +  \  ^ 


The  substitution  of  (2.13)  into  (2.12)  provides 


811 

1 

(2.14) 


1 

+  1 


*  (i  +  +  i 

/»<“  I  4'.  1 


r  )SS1 

iV 


4’,  TL)na*”l/xo,T1)^l  . 


SIV 


-6- 


Writing 


(2.15)  n(4  I  Mq,  T1)  *  S*^;  M,  Mq)  , 

and  using  the  reciprocity  principle  given  in  a  preceding 
paper  [ 1 0 ] ,  we  get 

(2.16)  S*  (r1 ;  /l,  Mq)  =  s*  (Ti  J 


and  (2.14)  becomes 


(2.17) 


<5 


U)(T1)X*(M,  T1)X*(M0,  Tx), 


where 

(2.18)  X*(M,  tx)  =  1  +  ^  M'  ^/T-  • 

Equation  (2.17)  is  equivalent  to  that  given  by 
Bellman  and  Kalaba  [2],  Busbridge  [6],  and  Ueno  [lO]. 

Hence  the  diffusely  reflected  intensity  emergent 
from  the  surface  r  =  t1  is  yielded  by 


(2.19) 


I<V  +  JO  =  !iZ  S*(Ti;  M,  Mq)  • 
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3.  TIME- DEPENDENT  CASE 

The  equation  of  transfer.  Consider  a  plane-parallel, 
inhomogeneous,  and  isotropically  scattering  atmosphere  of 
finite  geometrical  thickness  >  0  when  the  lower 

surface  z  =  is  illuminated  by  parallel  rays  of  time- 

dependent  intrinsic  flux  Fo(t),  where  F  is  a  constant 
and  0(t)  is  the  Dirac  delta  function  of  time  t. 

The  equation  of  transfer  appropriate  to  the  present 
case  is  written  in  the  form 

(3.D  M81  (z'azt;~ +  *(z)I  =  2a(z)  /j_  I(z’ 

where  1  and  cr  are  respectively  the  volume-attenuation 
and  scattering  coefficients;  c  is  the  speed  of  light. 

Equation  (3.1)  should  be  solved  subject  to  the 
boundary  and  initial  conditions 

(3.2)  1(0,  t,  +  (U)  =  0  and  I  (zx ,  t,  ~  M)  =  \  6  &  ~ 

(0  <  H  <  1,  0  <  t), 

(3.3)  I  (z,  0,  +  fi)  =  0  and  I(z,  0,  -  H)  -  0 

(0<M<1,  0  <  z  <  zi^  • 


The  equation  of  transfer  in  the  time-dependent  diffuse 


reflection  field  is  given  by 


1 

(3.4)  u||  +  i  |i  +  l(z)I  =  j<l(z>/  Kz,  t.  M')< iM' 

+  |t,(z)e-(zl-z)/Mo  5(t  -  ^  ), 

together  with  the  boundary  and  initial  conditions 

(3.5)  1(0,  t,  +  M)  =  0  and  I(z,  *»  _M)  =  0  (0  <M  <1,  t) > 

(3.6)  I(z,  0,  +  /l)  =  0  and  I(z,  0,  -ji)  =0  (0  <  /i  <  1 ,  O^z^Zj). 

The  emission  probability  distribution.  Let 
P*(M;  z,  zl5  t-s)  (0  <  M  <  1,  o  <  z  <  zv  0  <  s  <  t)  denote 
the  probability  that  a  photon  absorbed  at  the  level  z 
and  time  s  will  reappear  in  the  direction  M  at  time  t 
in  the  radiation  emerging  from  the  surface  z  =  z^ .  We 
assume  such  a  medium  of  dilute  gases  that  the  duration  of 
temporal  capture  of  photon  is  negligible  compared  with 
the  mean  free  time.  Furthermore,  it  should  be  mentioned 
that  the  emission  probability  is  considered  as  homogeneous 
with  respect  to  time,  because  of  the  stationary  optical 
properties  of  the  medium. 

The  expression  for  p*(M;  z,  z^ ;  t)  consists  of  two 
terms:  The  first  represents  a  term  of  direct  transmission 

without  scattering  in  a  time  interval  t,  and  the  second 
is  due  to  multiple  scattering  of  photons. 


Then  we  have 


z  z 

p  *(M;  z,  z1;  t)  =  a(z)e  ^zi  5  (t - -  ) 


(3.7) 


z-z  '  I  ^  -lz-z'1/fi'  dM-' 


Putting  z  =  z,  ,  we  obtain 


p *(m;  zx,  z1;  t)  =  a(z1)5(t) 


(3.8) 


2<Hl  ^dz'  ^p.toiz'.z^t  -  ^ 


The  scattering  function.  Let  the  truncated  generating 
function  of  the  emission  probability  distribution 
p z,  z^ ;  t)  be  denoted  by 

(3.9)  n(M  I  Mq)  zx,  t)  =  z,  z^ ;  t  -  )e“(zl"z)/^ 

Under  the  assumption  that  TI(M  I  HQ,  z^,  t)  has  the 
Markovian  property  of  the  conditional  diffuse  reflection 
probability  at  the  surface  z  =  z^,  an  integral  equation 
for  n(M  I  fJL0,  z1?  t)  is  directly  derived  from  the 
Chapman-Kolmogorof f  equation  in  somewhat  modified  form, 
namely 


dz  . 


-10- 


n(/x  I  flo,  ZjL  +  A,  t  +  £(i  +  jj-)) 


(3.10) 


=  ft  TT(/i'  I  u  ,  z.,  S)p(|U  1  jU'  ;  z  ,  t-S,  A)  dSd/i '  , 
0  0  ° 


where  A  lies  between  0  and  z^, 


jf(z1)A 

n(M  I  H0,  zlt  S)  =  nCi  I  40,  z ! j  S)(l - jj -  ) 


(3.11) 


+  A  •  p*(/i;  Zj,  z^ ;  s)  , 


and 


p(/i  I  do;  Z]L,  t-s,  A)  =  p*(/i;  z1?  zx;  t-s) 


(3.12) 


I  (z. ) A 


+  5(m  -  MD)5(t  -  s)(l - jj- 


). 


In  the  limit  as  A  -»  0,  we  get 


an 

8z. 


(/i  I  no,  zr  t)  +  (i  +  i-)[~  It  +  *(*!>]  n 


(3.13) 


=  p*  (M ; z^ ,  z^  ;  t)  +  "2  I  »z^ , s)p*  (M ;z^ , z  ;  t-s)dsjp 


Inserting  (3,8)  into  (3.13),  we  obtain 


BfWJWrfSS.-iT  iKrfalBSWPBSIHSBW  SWOT*?** 
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i-+  <ji  + jr’^lt +  <(z>|n  =  t'<z)lo(t> 

1  o 

(3.14)  +  i  /n(M  I  4’,  zx,  t)|4~  +  J  /W'  I  M0,  z,  t)^- 

+  4  ///' ^  (M  i  M '  ,z,s)n(M"lM  ,z,t-s)ds^-  Sr-  1  • 

4  0  0  0 


When  we  write 


(3.15)  II(M  I  M0,  z1,  t)  =  S(Z;L,  t;  |i,  /iQ) , 


(3.14)  becomes 


^  +  ‘<z»s 


(3 


.16)  =  jfejJltlt)  +|  /sCZj.tid.M'  )pr^  +  2  / SlZj^.t ;M",(i0)^S- 


4  ///< 


+  ^  j  J  J  S(z1>S;  M,M’) S(z1,t-s;  fi",M0) dS  |Jt-  1  • 


Equation  (3.16)  coincides,  in  the  case  of  isotropic 
scattering,  with  equation  (21)  of  [3] . 

The  requisite  diffusely  reflected  intensity  at  time 
t  is  provided  by 


I(z1(  t;  +  M)  =  |jj  S(z1,  t;  M,  MQ) 


(3.17) 
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The  principle  of  reciprocity.  In  a  manner  similar  to 
that  used  by  Biberman  and  Veklenko  [4],  we  shall  derive 
the  reciprocity  principle  as  follows: 

Let  f(x,  y)  be  the  joint  frequency  function  of  two 
one-dimensional  variables  £  and  t|,  which  is  expressed 
in  terms  of  the  corresponding  conditional  frequency 
function  of  f  (x  l  y)  ,  relative  to  the  hypothesis  T]  =  y, 
multiplied  by  the  marginal  frequency  function  f2(y). 

Then,  according  to  the  recipropity  relation  of  a 
random  process,  we  have 


(3.18) 


f  (x 


f (y  I  xjf-^x) 


We  write 


(3.19)  3(Zl,  t  -  s;  n,  HQ)  =  S(s,  Hq;  t,  M) , 

so  that  S (s ,  t,  #x)  (0  <  <  1 ,  P  <  <  1)  denotes 

the  probability  that  a  photon  absorbed  at  the  surface 
z  =  at  time  s  in  the  direction  will  be  found 

at  the  surface  z  =  z^  at  the  instant  t  in  the  direction 

M. 

Considering  that  x-state  and  y+state  variables  corres¬ 
pond  respectively  to  H  and  jit  variables,  we  get 


(3.20) 


S(^f  t;  H,  Mq)  =  S(zlf  t;  MQ,  H) 
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because  of  the  constancy  of  f^-function  (i  =  1,  2)  with 
respect  to  frequency  and  angular  arguments. 

In  this  context  it  should  be  noted  that,  when 
redistribution  of  radiation  energy  in  frequency  takes 
place  probabilistically  (i.e.,  noncoherent  scattering  is 
either  complete  or  real),  the  generalized  reciprocity 
principle  for  noncoherent  scattering  is  formulated  with 
the  aid  of  the  principle  of  detailed  balancing. 

By  use  of  (3,20),  (3.16)  reduces  to 


(3.21) 


=  a(z1)  ^  X(M,  Z] ,  s)X(mo,  z1 ,  t  -  s) ds , 


where 


(3.22) 


X(M, 


t)  =  6 (t)  + 


/ 


S (z, ,  t;  n,  M' ) 


d M’ 

JP 


When  the  optical  properties  of  the  medium  are  constant 
throughout  the  atmosphere,  (3.21)  becomes 


(3.23) 


+  ^  +  TroH^+l)S 


=  CU  ^  X(/Li,  T1,  u’)X(M0,  u 


-  u 


)  du ' 


where  T  =  z t,  t2  =  (l/c£)  ,  u  =  t/t2,  co  =  cr/1,  and 
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(3.24)  X(M,  T1,  u)  =  5  (u)  +  ^  /S(Ti»  u: 

In  the  above,  t ^  is  called  the  mean  free  time. 

In  a  later  paper,  with  the  aid  of  the  present 
probabilistic  approach,  the  generalized  principle  of 
reciprocity  for  noncoherent  scattering  will  be  treated 
Furthermore,  the  application  to  the  diffuse-reflection 
problem  of  radiation  in  spherical  and  cylindrical 
geometries  will  be  made. 
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